The Protogalactic Origin for Cosmic Magnetic Fields ^ 

Russell M. Kulsrud^'^, Renyue Cen^'^, Jeremiah P. Ostriker^'^, and Dongsu Ryu^'^ 

Received ; accepted 



^Princeton University Observatory, Princeton, NJ 08544 

^Dept. of Astronomy & Space Sci., Chungnam National University, Daejeon 305-764, 
Korea 

^e-mail: rkulsrud@astro.princeton.edu 
^e-mail: ryu@sirius. chungnam. ac.kr 
^e-mail: cen@astro.princeton.edu 
^e-mail: jpo@astro.princeton.edu 
^Submitted to the Astrophysical Journal 



- 2 - 



ABSTRACT 

It is demonstrated that strong magnetic fields are produced from a zero 
initial magnetic field during the pregalactic era, when galaxies are first forming. 
Their development proceeds in three phases. In the first phase, weak magnetic 
fields are created by the Biermann battery mechanism, acting in shocked parts 
of the intergalactic medium where caustics form and intersect. In the second 
phase, these weak magnetic fields are amplified to strong magnetic fields by 
the Kolmogoroff turbulence endemic to gravitational structure formation of 
galaxies. During this second phase, the magnetic fields reach saturation with the 
turbulent power, but they are coherent only on the scale of the smallest eddy. In 
the third phase, the magnetic field strength increases to equipartition with the 
turbulent energy, and the coherence length of the magnetic fields increases to the 
scale of the largest turbulent eddy, comparable to the scale of the entire galaxy. 
The resulting magnetic field represents a galactic magnetic field of primordial 
origin. No further dynamo action is necessary, after the galaxy forms, to explain 
the origin of magnetic fields. However, the magnetic field may be altered by 
dynamo action once the galaxy and the galactic disk have formed. 

It is first shown by direct numerical simulations, that thermoelectric 
currentsassociated with the Biermann battery, build the field up from zero to 
10~^^ G in the regions about to collapse into galaxies, by ^ ~ 3. For weak fields, 
in the absence of dissipation, the cyclotron frequency uJcyc = eB/mnc and 

+ x), where a; = V x v is the vorticity and x is the degree of ionization, 
satisfy the same equations, and initial conditions cucyc — cu — 0, so that, 
^cyclr,^) = ^(r, + x)- The vorticity grows rapidly after caustics (extreme 
nonlinearities) develop in the cosmic fluid. At this time it is shown that 
turbulence has developed into Kolmogoroff turbulence. Numerical simulations 
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do not yet have the resolution to demonstrate that, during the second phase, the 
magnetic fields are amplified by the dynamo action of the turbulence. Instead, 
an analytic theory of the turbulent amplification of magnetic fields is employed 
to explore this phase of the magnetic field development. Prom this theory, it is 
shown that the dynamo action of the protogalactic turbulence is able to amplify 
the magnetic fields, during the collapse time of the protogalaxy, by such a large 
factor that the power going into the magnetic field, must reach saturation with 
the turbulent power. There is as yet no analytic theory capable of describing 
the third phase. However, preliminary turbulence calculations currently in 
progress seem to confirm that the magnetic fields proceed to equipartition with 
the turbulent energy, and that the coherence length does increase to the largest 
scales. Simple physical arguments are presented that this is the case. Such an 
equipartition field is actually too strong to allow immediate collapse to a disk. 
Possible ways around this difficulty are discussed. 

Subject headings: cosmology: large-scale structure of universe - dark matter - 
hydrodynamics - magnetic fields - magnetohydrodynamics: MHD - methods: 
numerical - turbulence 
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1. INTRODUCTION 

How do cosmic magnetic fields arise? It has been the general belief that there is no 
natural mechanism to create a strong primordial magnetic field prior to the formation of 
galaxies. (For example, see Ruzmaikin et al. 1988.) For this reason, attention has been 
concentrated on the production of galactic fields after the galactic disks formed. It is 
assumed that some weak field, of order 10~^^ G, is initially present at the time of the 
beginning of the galactic disk, and this weak field is amplified to its present strength of a 
few /xG by dynamo action of interstellar hydrodynamic turbulence. (For a survey of the 
galactic dynamo theory see Ruzmaikin et al. 1988). But the theory for this amplification 
suffers from the problem of rapid amplification of small scale fields which inhibit the 
amplification of large scale fields (Kulsrud and Anderson 1992). In addition, the origin of 
the weak "seed field " is still open. 

Because of these difficulties it is of interest to reexamine the possibility that a strong 
primordial magnetic field can be produced in the pregalactic era. We have done that in this 
paper, and we come to the following picture of the magnetic field origin. 

Prior to the formation of galaxies, the universe is by no means uniform. Indeed, 
it is now known that small relative density perturbations, 5p/p, present at the time of 
recombination, grow to finite amplitude and form the present galaxies, clusters of galaxies, 
and other structures, at a variety of epochs. When 8p/ p ~ 0.1, shocks form and the resulting 
heated electrically conducting fiuid produces electric currents that generate magnetic fields 
even from field-free initial conditions. Such fields, which are quite weak (fa 10~^^ G), would, 
after compression into the galactic disk, produce the required seed field for the galactic 
dynamo. 

However, this is not the end of the story. The shocks also generate vorticity on 
all scales. This vorticity is strong enough that its vortex cells turn over and generate 
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Kolmogoroff turbulence. The energy spectra of this turbulence has the well known power 
law, and extends down to the viscous scale. (It is true that the medium is compressible. 
However, the shear motions in the turbulence are themselves incompressible. This is 
because the difference in frequencies between compressible and shear motions do not permit 
the sound waves to interact strongly with the incompressible shear motions. We can thus 
regard the turbulence as incompressible and of the nature considered by Kolmogoroff.) 

The eddies in the Kolmogoroff turbulence turn over at a rate that is fastest for the 
smallest eddy. This smallest eddy amplifies the magnetic field at a rate comparable to its 
turnover rate. The existence of a Kolmogoroff cascade from large eddies to small eddies is 
crucial because it enables the slowly rotating large eddies to drive rapidly rotating small 
eddies, that rapidly amplify the magnetic field. 

It will be shown that conditions are always such that during the collapse of a 
protogalaxy the smallest eddy turns over several hundred times. This would lead to an 
increase of the magnetic field energy by a corresponding number of powers of e, the base 
of natural logarithms, if saturation processes did not set in. Instead, the magnetic field 
saturates and comes into equipartition with the pregalactic turbulence. Such a saturated 
field is of sufficient strength to provide a primordial origin for galactic magnetic fields. 

If the dynamics of the magnetic field played no role, then the magnetic field would 
be very chaotic on small scales and the magnetic energy could be concentrated on scales 
smaller than the smallest eddy. However, as the field strengthens, the magnetic field 
on these smallest scales resists amplification, and the total magnetic field continues to 
strengthen only on scales comparable with that of the smallest hydrodynamic turbulent 
eddy. But later, when the field becomes even stronger, the smallest hydrodynamic eddies 
become suppressed due to the increasing drain of energy to the magnetic field. As a result, 
the spectrum of the turbulence becomes truncated at larger scales, and the scale of the 
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smallest hydrodynamic eddy increases. Eventually, only the largest eddies survive. During 
these later stages of saturation, only the magnetic energy at the largest scales is amplified, 
and the magnetic field eventually becomes coherent on these largest scales. 

Thus, during the pregalactic era the magnetic field goes through three phases. During 
the first phase, thermoelectric currents generated by shocks increase the magnetic field 
strength to a value of order 10~^^ G. This process is known as the Biermann battery 
mechanism after its discoverer (Biermann 1950). When this level of magnetic field strength 
is reached, the dynamo action of the turbulence becomes faster than the battery action, 
and the second phase commences. During this phase, the turbulence forms a Kolmogoroff 
spectrum down to the viscous scale. The smallest eddy at this scale does all the work in 
amplifying the magnetic field, and continues to amplify it until saturation sets in. Finally, 
as saturation is approached, the third phase takes over. During this third phase, the 
magnetic field approaches equipartition with the hydrodynamic turbulence, and becomes 
coherent on very large scales. 

To examine these three phases, we modified a standard code for the numerical 
simulation of the gravitational structure (Ryu et al. 1993). We first added the magnetic 
differential equation to the standard simulation equations to study the creation and 
evolution of the magnetic field, starting from a zero initial field. (Because the magnetic 
field is at first very weak, we neglected any magnetic forces in the numerical simulation, so 
the basic numerical solution for velocities and pressures was unaffected.) We found that 
the hydrodynamic turbulence does indeed create a magnetic field, and increases it to a 
field strength of about 10~^^ G, as expected for the first phase. However, the numerically 
simulated magnetic field did not increase beyond this point. The reason is that only a small 
part of the Kolmogoroff turbulence was captured in the numerical simulation. The smaller 
eddies, which would be the most important for dynamo amplification, were suppressed 
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by numerical resolution. The really important eddies (those on the smallest scale) were 
on a scale smaller than the grid size, and not seen at all. The only eddies captured in 
the numerical simulation turn over at too slow a rate to do more than overcome the the 
damping of the field by numerical resistivity. Thus, the simulation has too low a resolution 
to correctly represent the second phase of dynamo amplification. 

Therefore, the remaining conclusions concerning the enormous amplification of the 
magnetic field strength expected for the second phase can not be simulated until the 
numerical resolution is improved. To properly examine this amplification, we must employ 
an analytical method. To do this we assume that the Kolmogoroff spectrum inferred from 
the numerical simulation is real. This is plausible since the numerical evidence for the 
Kolmogoroff spectrum makes this assumption reasonable. The large scale eddies found in 
the simulations actually do turn over a sufficient number of times to produce a Kolmogoroff 
spectrum. Further, in the numerical simulations, where an approximately minus 5/3 power 
law was found over a limited range in wavenumber space, this range became larger when 
the numerical resolution was improved. Since the ideas of Kolmogoroff turbulence are now 
generally accepted, we felt secure in extrapolating the spectrum to the smallest viscous 
scales. 

Thus, assuming that the Kolmogoroff spectrum of turbulence is present during phase 
two, we employed a standard theory for the dynamo amplification of magnetic energy by 
hydrodynamic turbulence in the weak field hmit (Kulsrud and Anderson 1992). To evaluate 
the total amount of amplification we first scaled the total turbulent energy to the thermal 
energy. Prom the numerical simulations we know that these energies are comparable. Next, 
to find the scale of the smallest eddy we needed the plasma temperature and density. We 
normalized the temperature to the Jean's temperature. This gives a good approximation 
to the actual temperature during the time of the collapse of the protogalaxy. Using these 
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results we are able to make an estimate of the number of amplifications of the magnetic 
energy. We find that this estimate only depends only the normalization constants, which 
are of order unity, and on the the inverse square root of the galactic mass. For a baryonic 
mass in the range of a galactic mass, this estimate for the number of e-foldings of the 
magnetic energy is extremely large and demonstrates that the magnetic field will probably 
reach saturation with the turbulence. 

Our examination of the third phase, during which the magnetic field is strong and 
drives itself to coherence on a large scale, is still in progress (Chandran 1996b). However, 
physical arguments can be made which support the conclusion that the magnetic field 
should become coherent on a scale comparable with the size of the galaxy. 

Therefore, on the basis of the above considerations and calculations, we argue that the 
origin of galactic fields is primordial, i.e. predates the galaxies in their equilibrium state. 
Even if no magnetic field were present at the time of recombination, a magnetic field would 
be generated by the Biermann battery mechanism driven by the gravitational structure 
formation. Subsequently, during gravitational collapse of the protogalaxy to the galaxy, the 
residual turbulence will drive dynamo action that will increase the magnetic field strength 
to saturation. Finally, the strong magnetic field will smooth itself out on very large scales. 

The plan of the paper is as follows: In section 2, we describe the numerical simulation 
which shows the initial generation of the magnetic field, during the first phase, by the 
Biermann battery mechanism. In section 3, we apply the analytic theory of Kulsrud 
and Anderson (1992) to calculate the amount of dynamo amplification during the second 
phase. We show that this amplification is more than adequate to make the magnetic 
field sufficiently strong that the power into it reaches saturation with the hydrodynamic 
turbulence power. In section 4, we present the physical arguments that imply that during 
the third phase the magnetic field becomes coherent on large scales. In section 5, we present 
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our conclusions. 

2. THE NUMERICAL SIMULATION OF THE INITIAL GROWTH OF 

THE MAGNETIC FIELD 

The first phase of the evolution of the protogalactic magnetic field can be numerically 
simulated by modifying a cosmological hydrodynamic code normally used to simulate 
structure formation {e.g. Ryu et al. 1993). During this phase, the magnetic field is too 
weak to be of dynamic significance during the simulation. Thus, we include the magnetic 
field passively in the code because it does not affect the motions of the plasma. In the other 
words, we add the evolution equation for the magnetic field 



to the code, where v is obtained from the hydrodynamic part of the code without any 
magnetic forces. Thus, B can easily be followed without disturbing the rest of the 
simulation. 

The question arises concerning the initial value to take for the magnetic field in the 
modified simulation. According to equation (1), it must be initially nonzero if B is to be 
nonzero in time. A number of proposals have been made for the generation of a small initial 
magnetic field in the early universe. One of the more promising of these is the Biermann 
battery mechanism, which makes use of an extra pressure gradient term in Ohm's law 
(Biermann 1950). Employing this mechanism, we start with zero field at the beginning of 
the simulation, and include an extra term in the differential equation for the magnetic field. 
Thus, equation (1) is replaced by 





dB 

'dt 



V X (v X B) + 



(2) 



where Ug is the electron density, and Pe is the electron pressure. 
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If the ionization fraction, x, is taken constant in space, and the electron temperature is 
taken equal to the neutral temperature, then ne/nsi)- + x) = Pe/P- Here, ub is the baryonic 
number density. From this result, equation (2) reduces to 

f = Vx(vxB) + ^^^-^. (3) 

at e 1 + x 

where uih is the hydrogen mass, p = nsM/il + x) for a hydrogen gas. Equation (3) only 
involves the quantities which can be supplied from the hydrodynamic part of the code, so it 
is easy to incorporate the solution of equation (3) into the cosmological hydrodynamic code 
and to follow the evolution of B in time as the simulation proceeds. 

For any baratropic flow (p = p{p)), the last source term is zero because Vp is parallel 
to Vp. However, in general, for a real fluid in which curved shocks and photoheating can 
occur, Vp X Vp 7^ 0. 

Multiplying equation (3) by e/niHC, we get the equation for the cyclotron frequency 
ujcyc = eB/mnc 

ducvc ^ , N Vp X Vp 1 r7c „n , , , 

We have added a term to represent any resistive diffusion that may be present. 

Note that the ionization fraction enters into equation (3) through 1 + x, so even a very 
low ionization fraction is enough to generate magnetic fields. It is only necessary that there 
be enough electrons to carry the required current with a drift velocity relative to the ions 
vj:, less than their thermal velocity v^.^ i.e. B/A.t:L — n^evD/c < UeeVe/c where L is the scale 
size of variation of the magnetic field. This condition is satisfied by large margin in the 
numerical simulations. ( If this condition were not satisfied, an anomalous large resistance 
would develop due to plasma instabilities, and the magnetic field would be inhibited from 
growing.) 
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Equation (3) was incorporated into a numerical simulation for the large-scale structure 
formation in a standard cold dark matter (CDM) model universe with total Q = 1. The 
simulation was done in a periodic box with (32/i~^Mpc)^ volume using 128^ cells and 64^ 
particles from Zi — 20 to Zf — 0. The values of other parameters used are fib = 0.06, 
h — 1/2, and erg = 1.05. For the initial condition, we adopted the standard CDM power 
spectrum with n = 1, which is modified by the transfer function given by Bardeen et 
al. (1986). Thus, the simulation is basically the same as that reported in Kang et al. (1994), 
except a smaller box, and a smaller number of cells and particles was used. 

The actual equation for the magnetic field in comoving coordinates, which we solved is 

dB 1^ , ^, 1 ruHC 1 1 Vp x Vp 

— = -V X V X B - 2-B - —- ^ ^ ^ 5 

dt a ^ 'a Bote e 1 + , 

where Bg and to are the normalization constants for the magnetic field and time. This 

equation was solved by using the Total Variation Diminishing (TVD) scheme (Hartcn 1983). 

This is the same scheme as that employed for solution of the hydrodynamic equations (Ryu 

et al. 1993), so the accuracy was comparable. The constraint V • B = was enforced 

separately in every time step. We assumed x = 1. 

In Figure 1, the temporal evolution of the resulting magnetic field is plotted. The upper 
panel shows the volume-averaged (solid line) and mass-averaged (dotted line) magnetic 
energy density (B'^/Si:) as a function of z. The lower panel shows the volume-averaged 
(solid line) and mass-averaged (dotted line) magnetic field strength (B). Note that B ex h, 
and h — 1/2 was used. The magnetic field strength at first grows monotonically to the mass 
averaged value of order 10^^^ G by z ~ 2 — 3. In the core of clusters it reached 10^^'' G by 
z ~ 3. After this time, the value of the averaged field strength leveled off without further 
increase. It is believed that the saturation is due to the finite numerical resistivity inherent 
in the numerical scheme used to solve equation (5). 

The contours of the resulting baryonic density (p) and magnetic field strength (B) at 
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z — 1 are shown in Figure 2 and Figure 3. The shce shown has a thickness of 2/i~^Mpc (or 
8 cells). The upper panels show the whole region of 32 x 32/i~^Mpc, while the lower panels 
show a magnified region of 10 x 10/i~^Mpc. In Figure 2, the regions with density higher 
than the volume averaged value (0.06p) are contoured with contour levels (0.06p) x lO'^ and 
/c = 0, 0.1, 0.2, 2. Similarly in Figure 3, the regions with magnetic field strength 

higher than the volume averaged value (8 x lO^^^G) are contoured with contour levels 
(8 X 10~^^G) X 10*^ and A; = 0, 0.1, 0.2, . . . , 2. As expected, they are very well correlated. 
Since the magnetic field was mostly generated in the accretion shocks around the clusters, 
the high density core regions of the clusters have the strongest magnetic fields. 

It is interesting to compare the maximum rms value of S, or alternatively io^yc — lO^S, 
with the same rms mean for the vorticity, a; = V x v where v is the fiuid velocity. The 
maximum rms value for a; is ~ 10~^^s~^ around the clusters in the numerical simulation. 
This is equal to the cyclotron frequency of an ion in a magnetic field of ~ 10~^° G. In 
other words eB/niHC ~ a; in the clusters. This is not surprising since the equation for the 
evolution of —a; is identical to that for cUcyc = eB/mj^c, except for dissipative terms. 

By taking the curl of the equation of motion in the form 

— - V X (V X v) + -Vv^ = + i/VV (6) 

at 2 p 

where v is the kinematic viscosity, one gets 

— = V X V X u; - „ + uV^u. 7 
at p-' 

Now we see, on comparing equation (7) with equation (4), that if dissipative processes are 
ignored (conditions well satisfied except during the later stages of the simulation) , and if we 
assume that both cucyc and uj are initially zero, then we should have 
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a remarkable result. 

It must be appreciated that the Vp x Vp term is zero until some pressure is generated, 
since usually p is very small initially in the simulation. The generation of p generally 
happens in shocks where viscosity is certainly important. It can be argued that the jump in 
LUcyc and — <^/(l + x) across a shock should be equal since, if we could treat equation (7) as 
valid through the shock, the integral of ly'V^cu is probably small. Thus, a;cyc and cu satisfy 
essentially the same equation even in the shock. 

A check of the above relation is presented in Figure 4. The magnitudes of these two 
quantities are displayed on a logarithmic scale. Each point represents the two quantities in 
each cell. The magnitudes in one among eight neighboring cells were plotted. Here, h = 1/2 
was used again. If the relation in equation (8) holds exactly, all these points should lie on 
the line of unit slope. The deviation for small values is presumably due to the different 
dissipation rates which are not taken into account in the derivation of this relation. At 
larger values the correlation is much better, as is to be expected. The rough agreement of 
cUcyc and 00/(1 + x) at least for larger values tends to support the relation in equation (8). 

Eventually viscosity does become important and cu tends to saturate in mean square 
average. However, since the twisting of the magnetic field by the V x (v x B) term persists, 
one expects that B will continue grow. This fact is supported by Batchellor's discussion 
in his early paper (Batchellor 1950). It is thus indeed surprising that B seems to saturate 
at the same time with the same amplitude as cu does. Is it a coincidence that numerical 
resistivity becomes important at the same time that viscosity does? Since the effective 
viscosity is also numerical, this coincidence is not probably not so surprising. 
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3. THE ANALYTICAL CALCULATION OF THE EXPONENTIAL 
GROWTH OF THE PROTOGALACTIC MAGNETIC FIELD UP TO 

SATURATION 

Assuming the saturation of the magnetic field and uj are numerical, we can argue as 
follows with respect to what really think happens to the magnetic field: The value of the 
numerical viscosity is not much different from the physical viscosity, at least within a few 
orders of magnitude. But, the numerical resistivity is larger than the physical resistivity by 
great many orders of magnitude. Thus, it would be natural to assume that, in the absence 
of numerical resistivity and viscosity, the magnetic field should continue to grow, but uj 
should not. This supposition is strongly supported by Batchellor's argument (Batchellor 
1950). 

To numerically explore the beginnings of the second phase, in which the magnetic 
field strength grows exponentially, we plan to repeat the simulation under conditions in 
which the resistivity is smaller than the viscosity. The numerical resistivity and numerical 
viscosity are fixed by the grid size. However, if an additional viscous term is added to 
the simulation equations, with the viscosity coefficient clearly larger than the numerical 
viscosity, then this condition should be achieved. Such a simulation should enable us to 
confirm that the exponential growth of the magnetic field in phase two is real. 

The present numerical simulations have too low a resolution to follow the behavior of 
the magnetic field into the later phases of evolution, when dynamo action is important. 
Thus, we turn to the analytical theory of Kulsrud and Anderson (1992) to estimate the 
amount of the growth of the magnetic energy. For this purpose, we assume that during 
this time the turbulent motions are represented by Kolmogoroff isotropic homogeneous 
turbulence. The normalization constant in the Kolmogoroff spectrum will be determined by 
the numerical simulation. 
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Let us examine the three-dimensional power spectra found in the numerical simulations. 
The numerically determined spectra for (V • v)^, Pd(/c), and for (V x v)^, Pc{k), at ^ = 0, 
are plotted in Figure 5. At long wavelengths, the amplitude of the perturbations are small, 
so that hnear theory applies. In agreement with the results of Peebles(1993) and others, 
Pc{k) — > as A; — > 0, while Pd{k) follows the analytic theory expectation, Pd,{k) ~ . 
Thus, for very long wavenumbers, the (V x v)^ spectrum is very small, while that of the 
sonic turbulence is large. For these waves the shear turbulence has not yet gone nonlinear. 
(Since uj is driven by the Vp term as in equation (7), it cannot really develop until shocks 
form.) 

For wavelengths smaller than about 10h~^ Mpc, nonlinearities dominate, and we see 
that Pc{k) / Pd{k) ^ 1 . This spectrum is essentially the shear turbulence spectrum. The 
corresponding spectrum for sonic turbulence is also given on this plot. 

The shear turbulence peaks at a wavenumber of about one reciprocal Mpc . For values 
of k somewhat larger than this peak value, the spectrum follows a power law of k~^^^. 
(This power law behavior for the spectrum is expected for Kolmogoroff turbulence. This 
spectrum should develop from values of k for which the eddies have time to turn over once, 
since according to Kolmogoroff this turn over generates a cascade following the Kolmogoroff 
power law. Apparently the eddies at the peak have not yet turned over once, since the 
power law spectrum does not extend all the way to the peak. However, at slightly larger 
wavenumbers, the eddies have apparently turned over, as is indicated by the fact that the 
Kolmogoroff spectrum commences at these /c's.) The five thirds law does not continue to 
the smallest scales, but falls more rapidly with k as the grid scale is approached. This cutoff 
is expected for the numerical simulation because of the large numerical viscosity. 

To demonstratethe important result that that the real turbulent spectrum continues 
to smaller scales, one can compare the spectrum for a finer resolution simulation (see 
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Figure 5). It is found that the spectrum for the finer resolution extends to smaller scales, 
as is expected from the smaller effective viscosity if the turbulence is Kolmogoroff . . This 
is evidence that, in the real physical situation, there should exist a complete Kolmogoroff 
spectrum extending down to the much smaller true viscous scale. 

Although the numerical results are presented ai z = 0, the qualitative behavior should 
be the same as at earlier values of z. We apply the general Kolmogoroff theory to earher 
epochs when most galaxies were thought to have formed, and employ similar scalings for its 
amplitude. 

Let us introduce the one-dimensional spectrum for 

= y I{k)dk. (9) 
Then, if the turbulence is isotropic, we have 

(V X v)2 = 47ry P{k)k^dk = J k'^Idk, (10) 

so that AnP{k) = I{k). 

We assume that the Kolmogoroff spectrum starts at the peak, kmin, and extends down 
to some maximum k, k^ax, where it is truncated by physical viscosity. That is, we write the 
Kolmogoroff spectrum as 

^ - k^/3 ' ^^^) 
where Vq is the mean square turbulent velocity. 

Prom the analytic theory for the build up of magnetic energy given by Kulsrud and 
Anderson (1992) we have 

^ - 2;e.,. (12) 
where £m is the magnetic energy. 7 depends on the Kolmogoroff turbulence as 

2,./^*. (13) 
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where Aa;^ is the decorrelation rate at each k size. We take Aa;^ to be the eddy turnover 



rate kvk-, where Vk ~ \JkI{k) is the typical eddy velocity at the k scale. Combining these 
results we find 

27 = / ^/kIdk ^ voklll ^ ^ Ivokll^^kUL (14) 

if kmax ^ kmin- This is essentially the turnover rate of the shortest eddy with wavenumber 
k 

""max- 

Now, for a kinematic viscosity, u = Ivth, where / is the mean free path and Vth is the 
ion thermal velocity, we have 

"'mm \"'min'^ J \'-'th"'mm^ J 

Therefore, 

o„.2l3/2 ,1/2 

A more accurate value for 7 is derived in the Appendix. It is 

,1/2 3/2 

T= 1.236 (17) 

In deriving this formula we set k^in — ko, and assume the spectrum is zero for k < kg. 
The mean free path, Iq = IO^'^T^^/ub cm, and the thermal velocity Vi = 10^-\/7^ cm/sec. 
Tgy is the temperature in electron volts. The constants a and e are of order unity. They 
are introduced in the Appendix, and are related to the model for Kolmogoroff turbulence 
employed there. The quantity 27 is the exponential growth rate of the magnetic energy in 
a stationary fluid, and 7 is the growth rate of the root mean square fleld strength. 

In the expression for 7, we expect Vq to be related to the thermal energy ZusksT, 
where /c^ is the Boltzmann constant. Let 

^pvo = SPuBkBT, (18) 
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where (5 is the ratio of turbulent energy to thermal energy. (This ratio should be of order 
unity since behind a strong shock the ratio of kinetic energy to thermal energy is 15/16.) 
Thus, 

Let us imagine that the plasma is initially in the form of a sphere of radius R — 2'K/kQ 
and collapses to form a galaxy. (The collapse is more likely to be two-dimensional forming 
an oblate spheroid rather than a sphere. In this case R should refer to the shorter axis. 
However, the results on the amplification of the magnetic field should not be very different 
from those for a spherical collapse.) 

The magnetic field is also amplified by the collapse, as well as by the turbulence. We 
assume that the largest turbulent eddy has the same radius as the galaxy so — 27: /R. 
Then, if B is the rms value of the magnetic field we have 

We expect the sphere to collapse by a factor of order 2 in the dynamic time 

tn - 1 - X 10^^ .21) 
I'D — I = — -j= — W — s, yM) 

^jA-KGuBmHipD/pB) V PD 

where Pd/ Pb is the ratio of dark matter density to baryonic density. 

Thus, the amount of amplification during the collapse of the sphere by a single factor 
of two is of order 7^1). Substituting ko = 27i/R, and introducing the numerical values for 
lo, Vi and to, we get 

13'" 9.67 X lO'" 

Now, at the time of collapse, we expect that the temperature T should be comparable 
to the Jean's temperature. Therefore, we take 

ksRT^v = aCMniH, (23) 
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where M is the total mass of the sphere, dark matter plus baryonic matter. Let 
M = 2 X lO^Miip/j/psgm, where Mu is the baryonic mass of the galaxy in units of lO-*^^ 
solar masses, a is the ratio of T to the Jean's temperature, and is less than unity. During 
the collapse the temperature should remain isothermal because of cooling, or at least it 
should remain smaller than the Jeans temperature. For simplicity, we take a to be constant. 
Its actual effective value should be smaller than unity because of cooling. Thus, 



To form a more precise estimate of the factor by which BR'^ grows, we need the integral 



Consider as an example a uniform-density, isolated sphere of radius R. Let the initial 
radius be Rq, the initial value of 7 be 70, and the initial value oi to he tno- Let R — xRq. 
During the collapse of the sphere, the pressure of the baryonic matter is negligible so that 
the equation for the time dependence of x is 

d^x An ^ 1 

where po is the initial density. 

Integrating equation (26) we have 




(24) 




(25) 




(27) 



For simplicity we neglect this constant, so that 



dx 
'dt 




(28) 



From equation (24) we have that as time varies 7 ~ t 



-1 

D 



P 



,1/2 



~ X 



Using this 



scaling, we can now integrate 7 in time 




(29) 
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Substituting for ^yotoo from equation (24), we have 




(30) 



Taking Pb/pd 



1/20 and x 



1/2, we have 



109.7 (3^/^ 



(31) 



VMn(7eV2Q;i/2- 



The number of powers of ten by which R^B increases during the dynamic time, Aio, is 



4. THE EVOLUTION OF THE MAGNETIC FIELD AFTER SATURATION 

An increase in the rms magnetic field strength by this large number of powers of ten, 
clearly indicates that B should grow to reach saturation with the turbulent power. Just how 
the saturation proceeds is somewhat complicated, and deserves further analysis. However, 
from the results of Kulsrud and Anderson (1992), we can make some remarks as to how we 
think that the saturation should proceed. 

First, before saturation occurs, during the period when the magnetic energy grows 
exponentially, the wavenumber at which the magnetic energy is concentrated propagates 
to very small wavelengths. The peak wavenumber kpeak of the magnetic spectrum, Msik), 
increases exponentially, as kpeak ~ e'''*/^. The magnetic spectrum itself increases as k^l"^ up 
to kpeak-i and then falls off rapidly at larger wavenumbers. At any fixed k, the magnetic 
spectrum grows in time as e^'^^^^. The total integrated magnetic energy, 8m = J Msik), 
thus, increases as as e^'''* as it should. 



Second, the magnetic energy at scale k has an alfven frequency a;^(/c) = k^2SM/p 
(Chandran 1996a). As Sm increases and kpeak increases, there comes a time when the 
frequency ouAikpeak) exceeds the eddy turnover rate 7. When this happens, the eddy 




(32) 
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frequency become mismatched with the alfven frequency, and amphfication of magnetic 
energy at these scales is no longer possible. This means that for scales such that k > kcru, 
where 



the magnetic energy no longer grows exponentially. It is seen that as Sm increases, kcru 
decreases. It can be shown that, by the time that kcru becomes as small as k^ax, the 
wavenumber of the smallest turbulent eddy, the magnetic energy is concentrated at kmax- 

Third, at this latter time 2jSm is comparable to the turbulent power (ft! kpvl/2). 
After this time the magnetic energy stops growing exponentially and starts growing linearly 
in time. This transition happens because the drain on the turbulent power by the magnetic 
field is stronger than the drain by viscous dissipation. Hence, the hydro dynamic spectrum 
is truncated by the drain of energy into the magnetic energy, rather than by viscous 
dissipation. As a consequence kmax decreases at just such a rate that 7, which depends 
on kmax as k'^l^, decreases so as to keep always equal to the turbulent power. At 

the time when saturation commences, the cutoff wavenumber kcru — kmax- Subsequently, 
the wavenumber range in which the magnetic energy is predominant continues to be at 
kmax-i and this wavenumber becomes smaller and smaller. In the end, when essentially all 
the turbulent energy has been converted to magnetic energy, the magnetic energy is all 
concentrated near kmax ~ ^mm- 

To determine if this last phase is reached, we must ask whether, during phase two 
when exponential growth occurs, the magnetic energy at kmax reaches the value Esati given 



This is the value of magnetic energy at which the turbulent power into the magnetic field is 
comparable with the total turbulent power. (At this time, the turbulent power at smaller 
scales has stopped growing. The rate of growth of magnetic energy is (8/4)7, rather than 




(33) 



If -1 



(34) 



-22- 



27-) If this condition is satisfied, then we expect that all the turbulent energy will convert 
to magnetic energy. Thus, Ssat is a critical value for Em- At this magnitude and beyond, all 
the magnetic energy is in a band less than kmax- The energy in this limited range has grown 
only by the factor e^'^^^. Thus, during collapse of the galaxy, the important number is the 
number of e-foldings of the magnetic energy in the range ko to kmax- For total saturation 
it must be large enough for Sm to reach this critical value. Since number of e-folds of 
the energy in this limited range is (3/4) / jdt, the number of powers of ten by which the 
effective value of BR"^ increases during the collapse is reduced from that of equation (32), 
by the factor 3/8, that is to 

4/ _3, _ 17.87 P'/' 

~ 8'' ~ TK^wTv- ^^^^ 

If this number is large enough to amplify the magnetic field up to the saturated value such 
that 

^7^M = ^kminPVo, (36) 

then the remaining turbulence energy will be converted to magnetic energy at a linear rate, 
and the field will become coherent on the scale of the largest turbulent eddy, k^m- the 
other hand, if this number is too small to reach the saturated value of the magnetic field 
before the turbulence damps, then the resulting magnetic field strength is smaller, and the 
magnetic field is much more chaotic. 

For the determination of the critical number, A[q, we need values for a,P,a, and 
e. For a lower limit on A[q we may take a = 1, since we expect the temperature at the 
commencement of collapse to certainly be less than the Jean's temperature Tj. A realistic 
estimate of a, where aAcUk is the breadth frequency spectrum of the turbulence at fixed 
k, is about 2. (The definitions of these parameters are given just before equation (A4).) 
A value of e of 1 seems plausible where the parameter e is related to the rate of turbulent 
transfer (see equation (A14)). To find P we make use of the numerical simulation. P is the 
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ratio of turbulent kinetic energy to thermal energy. If we assume that the turbulent energy 
is one half of the total kinetic energy, then we can take the ratio of total kinetic energy 
density to thermal energy to be 2j3. This ratio depends on the ratio of the density to the 
mean density, p/ < p >. Its value taken from the numerical simulation, is plotted in Figure 
6. The appropriate value for p/ < p > at the time of collapse is 5 (Peebles 1993). Prom 
Pigure 6 we find from this value for the relative density that /3 ^ 5. Thus, substituting 
q; = 1,/? = 5,(7 = 2, and e = 1, into equation (35) we find that the number powers of ten 
by which B increases is 

A'lo = 30. (37) 

This number is certainly sufficiently large to show that B increase from the value 10~^^ G 
obtained from the Biermann battery mechanism, to the critical value of B for saturation 
with the turbulent power. (Por the value of 7 found from equation (24) with our choice of 
parameters, this critical field strength for B is about 10~^ G. If this critical value for the 
magnetic field strength is reached, then we expect equipartition with the total turbulence 
energy to result. This final value of B is about 10~^ G.) 



5. CONCLUSION 

In conclusion, we have described a mechanism for generating the observed cosmic 
magnetic fields that proceeds in three phases. 

In the first phase, the Biermann battery mechanism, driven by cosmic turbulence 
originating in shocks and caustic formation, creates a magnetic field from a zero initial 
field. This mechanism is found to produce a magnetic field whose strength is about 10~^^ 
G. In the second phase, Kolmogoroff turbulence amplifies this field by a factor sufficient 
to reach saturation with the turbulent power. (Neglecting saturation the factor of increase 
would be of order 10^°.) In the third phase, the magnetic field reaches equipartition with 
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the turbulent energy, yielding a field strength of order 10~^ G. During this phase it becomes 
coherent on the scale of the largest turbulent eddy. This scale is comparable with the scale 
of the entire galaxy. 

The final value of the magnetic field, which we predict to be in equipartition with the 
turbulent energy, is very strong. In fact, it is too strong to allow interstellar matter to 
flatten into a disk once it separates from the dark matter halo. The likely state of affairs 
at this time is that the sphere of interstellar plasma cannot immediately collapse to a disk. 
Probably, what happens is that star formation occurs removing a substantial fraction of 
the interstellar plasma from the magnetic field. The remnant of the plasma is then too 
light to hold the magnetic field, and the bulk of the remaining plasma is expelled back into 
the intergalactic or intra cluster medium carrying the field with it. The remaining field can 
then be compressed by the interstellar medium, and enriched by matter thrown off by the 
evolved halo stars. Such a process should leave some observational trace, which might be 
used to check this theory of the protogalactic origin of the cosmic field. 

It should be remarked that the dynamo action that we have examined has no a effect 
such as is usually postulated for mean field dynamo theory. These theories are kinematic 
and ignore any effect of magnetic forces. Thus, to produce a coherent field in the kinematic 
limit, some coherence in the turbulence is necessary. However, when the field becomes very 
strong, it generates its own coherence, since it is impossible to put a large energy into small 
scales. This is because the magnetic tension would become too strong if the magnetic field 
strength is increased on small scales. The magnetic field is perfectly capable of unwinding 
itself, and producing its own coherence, when it is sufficiently strong. The necessity for an 
a effect is purely a weak field result, only applicable in the kinematic limit. 

It is interesting that the triple phase process of creation of cosmic magnetic fields was 
suggested by Biermann himself many years ago(Biermann 1950). 
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FIGURE CAPTIONS 



Fig. 1 Temporal evolution of the magnetic field. The upper panel shows the volume- 
averaged (solid line) and mass-averaged (dotted line) magnetic energy density {B^/Stt) 
as a function of z. The lower panel shows the volume-averaged (solid line) and 
mass-averaged (dotted hue) magnetic field strength (B). 

Fig. 2 Density contours of a slice with a thickness of 2/i~^Mpc (or 8 cells) at z — 2. The 
contour lines with density higher than 0.06p are shown with levels 0.06 x 10^ and 
A; = 0, 0.1, 0.2, . . . , 2. The upper panel shows the whole region of 32 x 32/i~^Mpc, 
while the lower panel shows the magnified region of 10 x 10/i~^Mpc. 

Fig. 3 Magnetic field strength contours of a shce with a thickness of 2/i~^Mpc (or 8 cells) 
at z — 2. The contour fines with magnetic field strength higher than 8 x lO^^^G are 
shown with levels 8 10~^^ x lO'^ and k = 0, 0.1, 0.2, . . . , 2. The upper panel shows 
the whole region of 32 x 32/i~^Mpc, while the lower panel shows the magnified region 
of 10 X 10/i-^Mpc. 

Fig.4 Magnitude of a;/(l -|- x) plotted against that of LUcyc on a logarithmic scale. Each 
point represents the values in each cell. One among eight neighboring cells were 
plotted. The predicted relation is the forty-five degree straight line. The correlation 
is quite good for the larger values. 

Fig. 5 Three-dimensional power spectra for (V • v)^ (solid curves). 
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and for (V x v)^ (dashed curves). 
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produced by numerical simulations. The thick (solid and dashed) curves are from a 
power law (P^ = k~^) hydrodynamic simulation with boxsize L = 80/i~^Mpc and 512^ 
fluid elements, and the thin (solid and dashed) curves are from a power law (P^ = 
hydrodynamic simulation with boxsize L — 90/i~^Mpc and 288^ fluid elements. 
The L — 80/i~^Mpc simulation has a resolution half as that of the L — 90/i~^Mpc 
simulation. Comparison of the results from the two simulations indicates that the 
steep decline at the high k end of the power spectra is due to limited numerical 
resolution. Also shown is a straight dotted line which is parallel to the Kolmogoroff 
spectrum of logarithmic slope —5/3. 

Fig. 6 Ratio of the local kinetic energy to the thermal energy density in a numerical 
simulation, as a function of the local relative density p/ < p >. < p > is the mean 
density. If the turbulence kinetic energy density is one half of the total kinetic 
energy density, then this ratio is equal to 2/3, where /3 is the normalizing parameter 
introduced in the text as the ratio of turbulent energy density to thermal energy 
density. 
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A. A More Precise Value for the Growth Rate 7 



In this appendix, we evaluate the growth rate of the rms magnetic field for the 
Kolmogoroff spectrum of the turbulent energy per unit mass per unit wavenumber, I{k). 
I{k) is the one- dimensional spectrum defined in Eq. (9). We assume that / is zero for 
k < ko, and in the inertial range k <^ kmax where kmax is the cutoff wavenumber, I{k) is 
given by 



Here, Vq is the mean square value of the turbulent velocity. I{k) suffers a viscous cutoff at 
kmax- We make use of the formahsm of Kulsrud and Anderson(1992). 

In Kulsrud and Anderson, Vq/2 is given in terms of the k,u! spectrum, J{k,u!), 
the energy per wavenumber mode, per unit angular frequency, per unit mass in a box 
normalization of size L, 




(Al) 




(A2) 



while 7 is expressed in term of the quantity 




(A3) 



which involves the zero frequency harmonic of the turbulence. 



We assume the frequency spectrum at each k has a width (7(Aa;)fe, where 



(Aa;)fe = kvk, 



(A4) 



and (7 is a constant of order unity. Here, 




2k 



(A5) 



is the contribution from a range about k extending from k/2 to 2k. 




(A6) 
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so that 

{Auj)k = 1.1980''l^/\ (A7) 

Thus, we have 

J J{k, uj)duj = J{k, 0)a{Auj)k. (A8) 
Combining the above equation with Eqs. (A2) and (A3) and the definition of I{k), we get 

Uik) = , • (A9) 

Substituting Eq. (A7) in this equation, we get 

U{k) = 0.2086^-47^ (AlO) 

Now, in the kinematic hmit, according to Kulsrud and Anderson(1992), the total 
magnetic energy £m grows at the rate 

^ = 27^M, (All) 

where according to their Eq. (2.44) 

1 P^max _ „ 4-7r r^max , 

7 = - / kWk = — / k^Udk. (A12) 
Substituting Eq. (AlO) for U, we get 

7 = 0.8738 / ^ ' dk. (A13) 

Jko (7 

The KolmogorofT spectrum I{k) is given by Eq. (Al) for k not too near kmax- In order 
to find a more complete expression for / near k^ax, we carry out a dimensional analysis of 
the Kolmogoroff spectrum. 

We assume that, due to mode coupling, all the energy from a band k/2 to k is moved 
to the band k to 2k in a time l/e(Aa;)fe. Here, e is a constant of order unity. Thus, the fiux 
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of energy passing through /c, is 

$ = ^e{Auj)k f Idk = e(Aa;)fc|(22/3 - l)kl = 0.4406e(Au;)fcA;7. (A14) 

In the evaluation of the constant we have employed the power law approximation for the 
inertial range of I{k), Eq. (Al). 

Now, according to Braginski, energy is dissipated by viscosity at the rate 



\ dz ] ~^ \ dx 



+ 



(A15) 



where z is the prevailing direction of the rms magnetic field. The above result is correct as 
long as the gyroradius is small. Averaging this result over all directions, we get the energy 
dissipated per unit k is 



5 5 . 

Thus, equating the derivative of the flux $ to this viscous damping rate, we get the 
approximate equation for / 



0.4408^ [e{Au;)kkI] = 



(A17) 



The solution to this equation that agrees with Eq. (Al) is 



k 



4/3 



where 



,3/4,1/4 



kmax — 5.033- ^3/4 



Substituting Eq. (A18) into Eq. (A13), we get 

7 = Y^kl^^k^Lvo- 



(A18) 



(A19) 



(A20) 
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Further substituting Eq. (A19) for k^ax, we get 

3/2,1/2 

7 = 2.050^^ (A21) 

As given by Braginski, the ion viscosity can be expressed as 

77 = 2.75vilo, (A22) 
where Vi is the ion thermal velocity and and /q is the ion mean free path given by 

Vi = lO^Tey cm/s (A23) 

lo^lO'^Tl/uB cm. (A24) 

Here, Tg^ is the ion temperature expressed in electron volts and is the ion density. 
Substituting Eq. (A22) into (A21), we get 

3/2,1/2 

^ - 1-^36 ll'' (A25) 
aey\' Iq' _ 

which is Eq. (17). 
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